Modulation of intersubband IR absorption under intense THz irradiation 
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We analyze the modification of the intersubband absorption of electrons in quantum wells under 
intense THz irradiation. An expression for the induced current is obtained, based on the adiabatic 
approach and the resonant approximation. We predict the occurrence of a significant fine structure 
as well as the broadening and shift of the absorption under THz pump in a MW/cm 2 intensity 
range. 
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I. INTRODUCTION 

Strong transverse fields have long been known for mod- 
ifying confined states in quantum wells (QWs). The 
examination of the interband optical transitions under 
transverse fields, both static and high- frequency, is a con- 
venient method to study these modifications (see refer- 
ences in 0,0 and 0,3, respectively). The excitonic ef- 
fect and the modifications of both electron and hole states 
under transverse fields have to be taken into account for a 
quantitative description of the interband linear response, 
ft is also interesting to study the intersubband response 
under infrared (IR) excitation of electrons between the 
ground and the excited conduction band states which are 
placed in a transverse field. Whereas the electro-optic 
modulation of the intersubband transitions is well inves- 
tigated 5], the influence of an intense THz irradiation 
on such transitions is not investigated to the best of our 
knowledge. In this paper we treat theoretically the effect 
of the THz pump on the IR intersubband absorption. 

The confined electron states in a QW, of width d, 
subjected to a transverse electric field E u cosui, are de- 
scribed within the adiabatic approach, if ftw <C £21, where 
621 /h is the frequency of the intersubband transitions. 
Since the levels oscillate with a frequency w, the (n + 1 )- 
order intersubband transitions, with n THz photons and 
a single IR photon, take place resulting in a fine struc- 
ture of the absorption. At the same time, the shape of the 
absorption peaks is modified under the THz irradiation. 
In addition, one may consider the THz irradiation as a 
perturbation if \e\E U! d/2 < 621- Otherwise a numerical 
description of the electron states have to be applied. 

The calculations below are based on the one-particle 
density matrix equation linearized with respect to the IR 
field Eq exp(— iftt), while the THz irradiation is taken 
into account in the framework of the adiabatic approach. 
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FIG. 1: Band diagramms for symmetric (a) and non- 
symmetric (b) QWs under transverse THz irradiation. 
Dashed and dotted lines schematically show variation of levels 
and potentials under maximal THz field. 
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The broadening is described by the phenomenological ap- 
proach which takes into account the LO-phonon emission 
in the spectral region > ujlo, where lolo is the opti- 
cal phonon frequency. We have considered two cases: a 
symmetric rectangular QW, and a non-symmetric one. 
The last case may be realized by adding a transversal dc 
electric field E D to the symmetric case, as can be seen in 
Fig. 1(a) and 1(b), respectively. 



The paper is organized as follows. In Sec. II we derive 
the relative intersubband absorption under THz pump 
starting on the density matrix equation. The case of the 
perturbative approach is considered in Sec. Ill, while the 
results for the numerical description are discussed in Sec. 
IV. Concluding remarks and a list of the assumptions 
made are given in the last section. 
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II. INTERSUBBAND RESPONSE 

In this section, we obtain the averaged over the 
THz pump period relative absorption of the IR probe. 
The high-frequency addendum to the density matrix 
Apt(z, z') exp(— iSlt) is governed by the linearized equa- 
tion written in the coordinate-momentum representation: 

dAp pt (z,z') if - , 
^ h -^(hzt - h z >t - ?M)Ap pt (z, z ) 

+^(5h z -Sh z >)p pt {z,z') = (1) 

with the transverse coordinate z and 2D momentum p. 
Here the Hamiltonian and the perturbation operator, h z t 
and Sh z exp(— iSlt), are given by 



Pz 

2m 



eE^zcosut, 8h z 



SI 



-Env, 



(2) 



where v z — p z /m is the velocity operator and w z is the 
potential energy inside the QW. Within the hard-wall 
approximation one have to consider the interval \z\ < d/2 
with the zero boundary conditions. The density matrix 
under THz irradiation satisfies the periodicity condition: 

Ppt+27T/u(z,z') = p pt (z,z'). 

It is convenient to use the parametrically time- 
dependent wave functions, (p zt , determined by the eigen- 
statc problem: 



±d/2,t 



0. 



(3) 



which is defined in the interval \z\ < d/2. The expansion 
of the first and zero-order density matrices Ap pt (z, z') 
and p P t{z, z') over this basis gives us: 



Ap pt (z,z') = ^Ap pt (fc,fc')^t>l 

kk' 

p pt (z,z') = $> pt (M')¥>M (4) 



kk' 



and Eq. (1) in the fc-representation takes the form: 

dA Ppt (k,k') 



dt 



+ ~{e kt - e k > t - hSl)Ap pt (k, k') 



[&P P t(k,k')S k ' k - 5f A Ppt {k' ,k) 



£ [Shf Ppt (k' \k') - Ppt (k,k')5h k ' k ' 



(5) 



We have taken into account that S kk — —8 k k * . 

The induced current density is written through Ap as 
follows: 



e 

L 2 
2e 



,d/2 

2_~ \ dz lim (v zi + v z >)Ap pt (z 1 ,z' 1 ) 
p J -d/2 *u*i-^* 



= j2j2 v k'k(t)Ap pt {k,k') 



(7) 



pkk' 



where L 2 is the normalization area and Vkk' (t) = 
J-d^dztfiztVztPzt- For the case of resonant excitation 
between the first and the second levels, when hSl ~ e 2 ii 
Eq. (7) can be transformed into jt — wi2{t)Apt{2, 1). 
Here we have also performed the summation over 2D mo- 
mentum according to Ap t (2, 1) = {2/L 2 ) ^ p Ap pt (2, 1). 
Using Eq. (5) and supposing that only the ground level 
is populated, we obtain an equation for Ap t (2, 1) in the 
form: 

9Apt ^' 1} + i p2i(t) -SI- iv 2 i(t)} Ap* (2, 1) 

= lm V2lit)n2D > 



(8) 



where Sl2i(t) = {e^t — £it)/fi is the time-dependent inter- 
level frequency, 1^21 (t) = fu — S 22 + 5} 1 is the effective 
relaxation frequency with the phenomenological relax- 
ation frequency uq, and n2D is the electron concentra- 
tion. The stepped frequency uq increases in the active 
region, SI > lulo, due to the emission of LO-phonons 
with the frequency (jJlo- 

Introducing the time-dependent conductivity tensor 
Cfit = 3t/En and solving Eq. (8), we write: 



e 2 n 2 p 
hSl 



V2l{t) / dt'v 2 i(t') 



(9) 



x exp 



j — ij dr [0 2 i(r) — SI — if 2 i(r) 



The averaged over the THz period relative absorption is 
determined as follows p: 



dt 



4tt r 1 

£n = —/=Re / . 



Cf2t- 



(10) 



Thus, in order to calculate £q one needs to find f2 2 i (*) 
and ^21 (t) and to perform the integrations over time. 



III. PARABOLIC APPROACH 



where the non-adiabatic factor S kk and the perturbation 



Sh kk are given by: 



kk' 



5 kk 



d/2 



dz 



d/2 



dv k *t ,.fc- 
dt 



-<Pzt> 



Skf = %E a 



d/2 



d/2 



dz<p k Jv z <p k z 



(6) 



The perturbative solution of the eigenstate problem (3) 
can be written as follows: 



k k , e E u Zkh' k' . 



k'^k 



£k ~ £k' 



(11) 



, c , 1 \eE u Zkk'\ 2 
Ekt ^ £fc + eE u Zkk cosut + y j cos tut, 



k'^k 



£k - £k' 
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FIG. 2 

narrow and wide QWs. Solid line: HQ. 
(E = kV/cm). Dashed (dotted) line 
non-symmetric case (QWs under dc field E = 10 kV/cm) 



na 



(2) 



for 



where the matrix elements Zkk' and the energy level 
are determined from the zero-field eigenstate problem: 
(p 2 /2m + w z )<fz — £k<Pz, with the boundary conditions 
f\z\=d/2 — 0- The frequency O21 (t) in Eq. (9) is written 
as follows: 



^21 (t) = ^21 + cos uit + fl^i cos2wi, 



(1 2) 

where the coefficients il 2 i are obtained from Ekt given 
by Eq. (11). Note also, that Zkk = for the symmetric 
QW and the time-dependent contribution to (*) is 
cx cos 2ut for such a case. 

(2) 

The parabolic and linear dependencies of S7 21 and 

^21 on the applied THz field for both symmetric and 
asymmetric cases are plotted in Fig. 2. The calcula- 
tions are performed for narrow (d = 100 A) and wide 
(d = 250 A) Ga^4s-based QWs. In order to model non- 
symmetric QW we have introduced a transverse dc field 

E Q = 10 kV/cm besides the THz field E^cosujt. We 

(2) 

can see in Fig. 2 that the behavior of Q 2 i coefficients is 
exactly the same for the symmetric and non-symmetric 
cases in the region where the parabolic approximation is 
valid ( \e\E u d/2 < £21)- That is, for the narrow QW a 
parabolic bearing can be found up to very strong fields 
(E u ~ 300 kV/cm); whereas, for the wide QW, this ap- 
proximation is only valid up to E u ~ 20 kV/cm and 
beyond these regions Q2i(t) cannot be written by means 
of Eq. (12). For the non-symmetric case, Q 2 i coefficients 
show an almost linearly increase in the parabolic approxi- 
mation range for both narrow and wide QWs. After that, 
f2 2 i tend to a constant value. 

The non-adiabatic factors appearing in 1/21 (t) are given 



(12) 



-— y rzr sm 2uit = dfcfc sin 2ut. (13) 



Thus, V2i(t) can be written through Ais = 822 — <5ii- 
The velocity matrix element modifies weakly with time, 
V2i{t) — V21, because the time-dependent contributions 
(oc coswt, cos2Lat) are two order smaller than V2i- For 
the structures under consideration 1*21 — 4 x 10 7 cm/ s 
(d = 100 A) and w 2 i ~ 2 x 10 7 cm/s (d = 250 A). The 
non-adiabatic factor Av also shows a parabolic behavior 
which slightly depends on E a . For the region of parame- 
ters considered here, Ai//us varies from to 5 x 10 -3 and 
we have neglected these contributions below. 

Using the expansion of the exponential factors in (9) 
over the Bessel functions 6] and performing the integra- 
tions over time, we obtain £n for the symmetric QW in 
the form: 



47T e 2 v\v2i | 5 
■s/e he n 



-n 2 D 



E 



(fi 2 i -2kcu~n) 2 + . 



(14) 

In Fig. 3 we have plotted the spectral dependencies of 
£n, given by Eq. (14), under different THz pumps for 
the narrow symmetric QW. The calculations are per- 
formed for a multiple quantum well structure with ten 
decoupled 100 A wide GaAs QWs and for THz quanta 
energy values hoj = 0.5 meV and 1.5 meV. We have 
used two broadening energy values hv ~ 0.66 meV (cor- 
responding to a relaxation frequency v — 1 ps _1 ) and 
hv = 1.5 meV. One can see that, for fields well below 
the parabolic approximation limit, the relative absorp- 
tion shows a monotonous decreasing while spreads to 
higher frequencies. For high fields, beyond 70 kV/cm, 

a clear multi-peak structure appears because the strong- 

(2) 

modulation case, when f2 21 /4cj > 1, can be realized in 
the framework of the perturbation approach. The rela- 
tive height of the different peaks for this structure de- 
pends on the quanta energy value, as can be seen from 
Fig. 3, due to the u> contribution to the Bessel functions 
and the denominators in Eq. (14). Panels (a) and (c), 
with 0.5 meV, show higher lateral peaks than the central 
one, whereas panels (b) and (d), with 1.5 meV, have a 
clear maximum in a central peak. The relaxation energy 
not only affect the broadening of the peaks but also the 
relative absorption and so, panels (c) and (d) with 1.5 
meV show wider and lower peaks than panels (a) and 
(b) in which the fine structure is more evident and the 
relative absorption is about twice bigger than that of the 
other cases. 

There are more sums over k for the non-symmetric QW 
due to the additional term with coefficient fi^i • I n this 
case the parabolic approximation leads to 



4-7T e 2 v\v2i\ 2 
\fi he fl 



-n 2 D 



E 

kki A/c— — 00 



(15) 
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FIG. 3: Spectral dependencies of £n versus hQ. for the sym- 
metric case under THz field strengths: = kV/cm (solid 
line), E u = 50 kV/cm (dashed line), E w = 70 kV/cm (dot- 
ted line) and = 100 kV/cm (dot-dashed line). Panel (a): 
hui — 0.5 meV, hv — 0.66 meV; panel (b): huj — 1.5 meV, 
hv — 0.66 meV; panel (c): hu) = 0.5 meV, hv — 1.5 meV; 
panel (d): fto; = 1.5 meV, hv — 1.5 meV. 



FIG. 4: Relative absorption for non-symmetric 100 A QW 
case under E — 10 kV/cm and = (solid line), = 50 
kV/cm (dashed line), E^ = 70 (dotted line) and E w = 100 
kV/cm (dot-dashed line), faj = 1.5 meV and hv = 0.66 meV. 



IV. NUMERICAL DESCRIPTION 



^ (¥) ^~2A fe ^+Afc ■/* (¥) 

[f7 2 i-(2fc-fci)w-f7] 2 + z/ 2 



Fig. 4 shows the relative absorption of non-symmetric 
QW for E Q = 10 kV/cm and different E u values. In or- 
der to obtain a higher relative absorption and a clearer 
peak structure, we have used parameters corresponding 
to panel (b) of Fig. 3 (hut = 1.5 meV and hv = 0.66 
meV). The behavior of the relative absorption is similar 
to that of the corresponding to the symmetric case for low 
fields, but the relative absorption decreases faster and, for 
fields beyond 70 kV/cm, when the multi-peak fine struc- 
ture becomes evident, there are clear differences between 
this case and the symmetric one. Looking at the 100 
kV/cm case in panel (b) of Fig. 3, we can see the symme- 
try of the relative absorption multi-peak structure with a 
central maximum and two lateral satellites. On the other 
hand, Fig. 4 shows a clear non-symmetric fine structure, 
with a pronounced central maximum located at the same 
energy than the corresponding to the symmetric case. 
However, lateral satellites show a strong asymmetry. At 
first glance it seems that a new plateau appears in the 
low energy side. Actually, for E u = 100 kV/cm, both 
symmetric and non-symmetric cases present seven rela- 
tive maxima as can be proved by diminishing the relax- 
ation broadening beyond realistic values. In other words, 
the absorption multi-peak structure is masked in part by 
both the THz quanta and the effective relaxation energy 
values. 



As we have already noted, the perturbation approach 
is valid under the condition \e\E u d/2 < e 2 i- If the THz 
pump is stronger, one needs to perform a numerical so- 
lution for the cigenstate problem (3) and to integrate 
Eqs. (9, 10) numerically taking into account the time- 
dependent velocity matrix element. Such a case can be 
realized in a wide QW without the requirement of us- 
ing high fields. Therefore, we will consider below 250 A 
wide GaAs-based QWs. In order to check the validity of 
the parabolic approximation, Fig. 5 shows the relative 
absorption calculated both by means of this approach 
and numerically. Calculations are performed for a mul- 
tiple QW structure with ten decoupled QWs and for a 
THz quanta energy fku — 1.5 meV. Due to the low val- 
ues of the energy range in this case (around the optical 
phonon energy), we consider a step- like function for the 
relaxation energy taking a smaller value hv = 0.5 meV 
in the passive region, when Ml < 35 meV, and a bigger 
one hv = 1.5 meV for higher energy values. One can 
see a very good correspondence between the two meth- 
ods for fields up to 30 kV/cm while, for higher fields, the 
parabolic approach clearly fails. 

The results obtained by means of the numerical de- 
scription for spectral dependencies beyond E w = 30 
kV/cm are represented in Figs. 6 (symmetric QW) and 
7 (non-symmetric QW). As in Fig. 5, calculations have 
been made for ten QWs and for the same THz quanta 
energy and step-like relaxation energy. For higher THz 
fields the asymmetry appears (even for the symmetric 
QW). Together with the decreasing and spreading of the 
relative absorption, the initial peak for E u = kV/cm 
gradually splits in several peaks as the field intensity in- 
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FIG. 5: Relative absorption £n obtained through numerical 
calculations (solid line: — 30 kV/cm, dotted line: E u = 50 
kV/cm) and by means of the parabolic approximation (dashed 
line: E u = 30 kV/cm, dash-dotted line: E u = 50 kV/cm) for 
the symmetric case and wide QW. 
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FIG. 7: Non-symmetric QW spectral dependencies of £n ver- 
sus fitl for Eo = 10 kV/cm and E^ = 30 kV/cm (solid line), 
50 kV/cm (dashed line), 70 kV/cm (dotted line). 



V. CONCLUSIONS 
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FIG. 6: Symmetric wide QW spectral dependencies of £n 
versus hQ, for E^ = 30 kV/cm (solid line), 50 kV/cm (dashed 
line), and 70 kV/cm (dotted line). 



creases, the number of peaks depending on the inten- 
sity. The non-symmetric QW well case shows two differ- 
ent regions not only due to the stepped relaxation but 
also because the field runs twice per cycle the region be- 
tween and E Q — E u and only once the remaining region. 
Thus, the multi-peak fine structure is partially hidden for 
hfl > 35 meV due, once again, to the higher effective re- 
laxation energy. 



We have studied the modifications of the intersubband 
absorption in a multi-QW structure caused by a strong 
THz irradiation. Experimentally, such an irradiation can 
be achieved by using free-electron or gas lasers with an 
energy density in the MW/cm 2 range (see applications of 
these lasers to study a heterostructure response in || or 
0, respectively). Results show a significant fine struc- 
ture of the absorption peak due to the (n-(-l)-order inter- 
subband transitions, with n THz photons and a single IR 
photon. In addition, a strong modification of absorption, 
which consists on a noticeable broadening of the zero- 
field peak and a shift towards higher energy values, is also 
demonstrated. Since the relative absorption calculated is 
weak enough, one can measure photoconductivity. 

Next, we discuss the assumptions used in the above 
calculations. Within the hard-wall scheme we have disre- 
garded the underbarrier penetration. In general, for deep 
and wide decoupled (or weakly coupled) quantum wells, 
underbarrier penetration slightly modifies wave functions 
and energy levels position. In the present case and due 
to the relatively big interlevel distance (£21 ~ 140 meV 
for narrow QW and £21 ~ 30 meV for wide QW) lit- 
tle changes in the level positions do not affect essentially 
results. For the same reasons we have not taken into ac- 
count possible contributions of the Coulomb renormal- 
ization neglecting the depolarization and exchange ef- 
fects. We have also used a phenomenological homoge- 
neous broadening energy hi> ~ 0.5 — 1.5 meV. These 
assumptions are generally accepted and a possible im- 
provement will not change essentially the obtained re- 
sults. We have also neglected the interlevel redistribu- 
tion under THz pump because the interlevel energy £21 
is bigger than 20 meV while the THz quanta energy hw is 
around 1 meV. One can neglect a THz field effect on the 
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damping because the influence of such a field on the wave 
functions is not very strong. Thus, a phenomenological 
broadening should not be essentially dependent on the 
THz pump. 

To conclude, the present work has shown the possibil- 
ity that an essential modification of the IR intersubband 



response takes place when an intense THz irradiation is 
applied. We expect that the results obtained encourages 
researches to carry out experiments in this direction. 
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